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The Dyson-Schwinger equation for the ghost-gluon vertex of the Hamiltonian approach to Yang- 
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gluon vertex is IR finite but IR enhanced compared to the bare one by 15% to 25%, depending on 
the kinematical momentum regime. 
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I. INTRODUCTION 

In recent years there have been intensive studies of continuum Yang-Mills theory in the non-perturbative regime. 
Most of these investigations were carried out either in Landau gauge, using Dyson-Schwinger equations (DSEs) 
(for reviews see Refs. [3-|3|) or functional renormalization group (FRG) flow equations (for reviews see Refs. 
|(|), or in Coulomb gauge, using DSEs 0j- In addition, variational [si— Tl5j| or FRG methods [n| were used in the 
Hamiltonian formulation of Yang-Mills theory in Coulomb gauge. All those various approaches have in common that 
the infrared (IR) sector of the theory is dominated by the ghost degrees of freedom, which has been referred to as 
"ghost dominance" . For this reason the ghost-gluon vertex is of crucial importance in these approaches. Inspired by 
Taylor's work [l7[ showing the non-renormalization of the ghost-gluon vertex in Landau gauge, in all the approaches 
mentioned above it was tacitly assumed that the ghost-gluon vertex is bare. In Ref. [l8| a semi-perturbative calculation 
, of the ghost-gluon vertex was carried out and it was found that its dressing is small even in the IR, a result which 
is also found on the lattice 0, |2(j. Although Taylor's proof of non-renormalization formally applies also to the 
ghost-gluon vertex in the Hamiltonian formulation of Yang-Mils theory in Coulomb gauge, the results for the dres sing 
of the ghost-gluon vertex obtained in the functional integral formulation of Yang-Mills theory in Landau gauge [lq 
cannot be assumed to remain valid also in the Hamiltonian approach in Coulomb gauge. Furthermore, recently it 
'— — was found within the FRG-approach in Landau gauge that the dressing of the ghost-gluon vertex becomes crucial at 
high temperatures [2l|. Since the high-temperature limit of the 4-dimensional Yang-Mills theory is essentially the 
3-dimensional Euclidean Yang-Mills theory and the latter provides an approximation to the Yang-Mills vacuum wave 
functional in 3 + 1 dimensionj^l(22j , the dressing of the ghost-gluon vertex of the Hamiltonian approach in Coulomb 
gauge should be expected to be also substantial. Therefore, in the present paper we investigate the ghost-gluon vertex 
of the Hamiltonian approach to Yang-Mills theory in Coulomb gauge. We will solve the DSE for the ghost-gluon vertex 
using as input the non-perturbative ghost and gluon propagators previously obtained in the variational approach. 

In Sec. [IT] we briefly summarize the basic ingredients of the Hamiltonian approach to Yang-Mills theory in Coulomb 
gauge and also present the results obtained for the ghost and gluon propagators. In Sec. lIIII we give a short derivation 
of the DSE for the ghost-gluon vertex, and introduce the truncation scheme. Our numerical results are presented in 
Sec. IIVI Some concluding remarks are given in Sec. [V] 
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^ ■ II. EQUAL-TIME PROPAGATORS IN COULOMB GAUGE 

U 

The Hamiltonian approach to Yang-Mills theory in Coulomb gauge is based on canonical quantization in Weyl gauge 
A$ — 0, and results in a Schrodinger equation, which has to be solved for the vacuum wave functional ip[A] = (A\0) 
of the transverse gauge field diAf — 0. Once i/)[A] is known all static (time- independent) Green's functions can be 
evaluated. In Refs. [T3, EH the Yang-Mills Schrodinger equation was solved in an approximate fashion using the 
variational principle and assuming Gaussian-type trial wave functionals. 



1 More precisely, the functional integral of 3-dimensional Euclidean Yang-Mills theory in Landau gauge can be interpreted as the functional 
integral of the vacuum expectation value of the Hamiltonian approach in Coulomb gauge in 3 + 1 dimensions with a vacuum wave 
functional given by tp[A] ~ exp(— ^ Sym[A\) , where S Y m[j4] is the classical action of 3-dimensional Euclidean Yang-Mills theory. This 
wave functional was shown to provide a decent approximation to the true Yang-Mills vacuum wave functional in the mid-momentum 
regime [23ll , 
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In the Hamiltonian formulation of Yang-Mills theory, Coulomb gauge can be implemented in the expectation value 
of any functional K[A] of the (spatial components of the) gauge field A by the Faddeev-Popov method, which results 
in 

(K[A})= f VAJ A \^[A]\ 2 K[A]. (1) 

Here, Ja — Det(G^ 4 1 ) is the Faddeev-Popov determinant with 

G A lab (x, y) = (~S ab 8 2 - gr h A1{*)d t )5(K - y) (2) 

being the Faddeev-Popov operator. Furthermore, g is the coupling constant and f acb are the structure constants of 
the su(N c ) algebra. The functional integration in Eq. (fTJ) extends over transverse field configurations restricted to the 
first Gribov region f2 or, more precisely, to the fundamental modular region [24]. Moreover, we assume that the wave 
functional ip[A] is properly normalized, such that (1) = 1. 

To simplify the bookkeeping we will use the compact notation ^^(xi) = A(l) and assume that a repeated label 
means summation over colour and spatial indices along with integration over the spatial coordinates, 

A - B = A{l)B{\) = J d 3 x A°WS?(x). (3) 

We use the same convention for indices referring to the ghost field except that the label "1" represents only the colour 
index a\ and the spatial coordinate xi. 

The gluon propagator D and the ghost propagator G are defined by the expectation values 

D(l, 2) = (A(1)A(2)) , G(l, 2) = (G A (1, 2)) . (4) 

In momentum space we express the gluon propagator by the gluon energy fl(p) 

<4 Q (p)4(q)) =: ^(2^(p + q)^|, (5) 

where Uj(p) = Sij — PiPj/p 2 is the transverse projector in momentum space. Furthermore the ghost propagator can 
be represented as 

G-*(p) = 5 ab ^, (6) 

where d(p) is the ghost form factor. Assuming the so-called horizon condition (i~ 1 (0) = 0, the results obtained with 
Gaussian-type wave functionals [13] show an IR diverging gluon energy fl(p) which can be fitted with Gribov's formula 



W = V p2 + "pf' (7) 

with a mass parameter m 2 A ~ 0.6 a c (for N c — 2), where a c is the Coulomb string tension, i.e. the coefficient of the 
linear term of the non-Abelian Coulomb potential. These results compare favourably with recent lattice calculations 
[26j : in particular, the infrared regime of the gluon energy is correctly reproduced. 



III. DERIVATION OF THE GHOST-GLUON VERTEX DSE 



In the functional integral formulation of Yang-Mills theory in Landau gauge the DSE for the ghost-gluon vertex 
has been known for quite some time, see e.g. Ref. [27j |. In the Hamiltonian approach in Coulomb gauge the DSE for 
the ghost-gluon vertex was derived in Ref. |15| . We briefly summarize this derivation in order to fix our notation. 

The Faddeev-Popov operator Eq. © can be inverted to give the following operator identity for the ghost operator 
G A 



G A (1, 2) = G (l, 2) - G A (1, 4)A(5) r (5; 4, 6) G (6, 2). 



(8) 
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FIG. 1. Diagrammatic representation of the DSE (|f 3[) . Small filled dots represent connected Green's functions, and small 
empty dots proper vertex functions. 



Here the bare ghost-gluon vertex Tq is defined by 



5G- A l {2, 3) 



rod; 2, 3) = -^, (9) 

and agrees with the lowest order term in the perturbative expansion of the full vertex T defined by 

(A(3)G A (1, 2)> =: -15(3, 3') G(l, V) G(2' , 2) f(3'; 1', 2'). (10) 

Taking the v.e.v. of Eq. ((5J) and using Eq. (flU)) one obtains the DSE for the ghost propagator G, which reads 

G-\l,2) = G \l,2) - f (3; 1, 4)£»(3, 3')G(4, 4')f „(3'; 4', 2). (11) 

If Eq. © is multiplied by a gauge field A(3) before taking the expectation value, we are led to 

(G A (l, 2)A(3)) = -f (5; 4, 6) G (6, 2) (G A (1, 4)A(5)A(3)) , (12) 

where we have used (A) — 0. For the left-hand side of Eq. (]12[1 the relation (ITOl) can be used, while the expectation 
value on the right-hand side can be expanded in terms of vertex functions and propagators [15j. By using the DSE 
([TT]) for the ghost propagator, the DSE (fT2"j) for the ghost-gluon vertex becomes 

f (1; 2, 3) = f (1; 2, 3) + f (1; 4, 5)G(4', 4)G(5, 5')f (6'; 5', 3)£>(6, 6')f (6; 2, 4') 

+ T(l, 4, 5)£>(4, 4')r>(5, 5')f (4'; 2, 6)G(6, 6')f o(5'; 6', 3) (13) 
- f (1, 5; 2, 4)Z>(5, 5')G(4, 4')f (5'; 4', 3), 

where r(l,4, 5) is the three-gluon vertex, and r(l,5;2, 4) is the ghost-gluon scattering kernel. Equation ([13]) is 
represented diagrammatically in Fig. [TJ Note that the vacuum wave functional does not explicitly enter this equation, 
but only implicitly via the various propagators, in particular the gluon propagator. 

Equation (1131) must be truncated to be feasible; in a first step we will discard the proper four-point function. Even 
then, a fully self-consistent solution of the resulting equation together with the DSEs for the propagators is still very 
expensive and beyond the scope of this work. To get a first estimate of the size of the dressing of the ghost-gluon 
vertex we will keep the full, non-perturbative propagators but approximate the vertices in the loop terms by their 
bare form. If the usual assumption of a bare ghost-gluon vertex is justified, the corrections to the bare vertex from 
the one-loop terms should turn out to be small. 

The bare three-gluon vertex is given by the three-gluon kernel 73 in the exponent of the vacuum wave functional 

m 

ip[A] = cxp \-\ujA 2 - -^IzA 3 — . • . 1 (14) 



2 3!2 

which is found from the variational calculation to be given by [l5j 

b 2/°*T„ fc (p >q> k) 

7 « fe(p ' q,k) -fi(p) + n( q ) + n(k)' (15) 

where T is the Lorentz structure of the three-gluon coupling in the Hamiltonian 

Tijfe(p,q,k) = ig[Sij(p — q)k + $jk{q - fe)< +Ski(k-p)j]. (16) 

Note that 73 Eq. (TT5|) is the bare part of the three-gluon vertex T3 in the sense that it is the leading term in the DSE 
for T3. However, 73 Eq. (|15p is not the lowest-order perturbative vertex, which is obtained from Eq. (|15l) when the 
gluon energy Sl(p) is replaced by its perturbative counterpart |p|. 
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After implementing the truncations explained above and extracting the colour structure, Eq. (I13|) reads in momen- 
tum space 



~ N f 

I\(k; p, q) = i 9 %(k) Pi - %(k)-^ J &£ (igp m ) G{£) (iglj) G(£ + k) (ig(£ + k) n ) 



20(-€-p) , , 

(17) 

, ,,^c /^w. ^ttmW 2T 3 - mn (M,-l-k) t nfc (l + k) ,., /j\ \ / p v V ^ 

where we have introduced the notation d£ = d 3 f/(27r) 3 . In Eq. (|T7)) k, p, and q are respectively the momenta of the 
gluon and the incoming and outgoing ghost, and momentum conservation k + p + q = 0is implicitly understood. 
We parameterize the full ghost-gluon vertex by a dressing function h as 

r,(k;p,q) = igt ij Q S .)p j [l + h(k 2 ;p 2 ,q 2 )]. (18) 

Let us stress that in the present Hamiltonian approach the Coulomb gauge condition is exactly implemented, and 
the functional integral of the scalar product of the Hilbert space is strictly restricted to (spatially) transverse gauge 
fields. Therefore, in the present case the ghost-gluon vertex cannot develop a longitudinal part. 
Contracting Eqs. (fTT)) and ([T5]) with pi and dividing both sides by PiPjUjO*) we obtain^ 

h(k 2 ;p 2 ,q 2 )=h(k 2 ;p 2 ,q 2 )+gl 2 (k 2 ;p 2 ,q 2 ), (19) 

where 

V[i-(k- P ) 2 ]7 JK ' J t 2 (i + k) 2 n(£-p) v ' 

is the contribution of the diagram with three ghost-gluon vertices, and 

N c f S£ P l p ] (p + k) k t u (k)t ]m (£)t kn (£ + k)T lmn (k,£, -I - k) d{£ - p) (2Qb) 



V[i-(k- P ) 2 ]7 n(€)fi(€ + k)[n(k) + n(£) + n(£ + k)] (£-p) 

is the contribution of the diagram containing the three-gluon vertex; the latter is multiplied by the coupling constant 
g due to the parameterization of the ghost propagator, see Eq. ©. It should be remarked here that the integrals 
Eqs. (|20[) are UV finite and need not to be renormalized. In the subsequent numerical calculations we will consider 
two classes of kinematic configurations, namely 

Mp 2 ;pW) and h(xp 2 ;p 2 ,p 2 ), (21) 

where x is restricted to the interval x G [0,4] due to momentum conservation. 

Before concluding this section we investigate the IR limit of the ghost-gluon vertex. In the IR the ghost and gluon 
propagator can be parameterized as 

rnft m\ +a 

with f3 > to fulfil the horizon condition d _1 (0) = 0. These IR exponents obey the sum rule 1 + a = 2/3 [l2l l28l l29j. 
Inserting the IR ansatzes Eq. (|22|) into Eqs. (|20l) we find that in the limit of vanishing momenta the form factor of 
the ghost gluon vertex approaches the finite expression 



fe(0;0,0)= — \ 2 — ■ (23) 
24(1 + p)-K z \m A J 

A fit to the numerical data [13( for the propagators with N c — 2 and the solution (3 — 1 yields m c ~ 4.97-^0^ while 
niA was given below Eq. (J7). Plugging these values into Eq. (f2"3f yields an IR value of 0.174 for the ghost-gluon- vertex 
form factor h. The numerical data shown in the next section confirm this value. 



2 The limit p ■ k — ► ±1 is finite. 
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IV. NUMERICAL RESULTS 



For the ghost and gluon propagators we use the results obtained previously in the variational approach [13( as 
input. In particular, for the gluon pr opagato r we use the Gribov formula [Eq. ([?])]. while for the ghost form factor 
[Eq. we use the parameterization [15|, |3(| 

d(pH V^ + ^^)' (24) 

with a = 4.97 and c = 16. The factor g has been absorbed in the ghost propagator as explained in Rcfs. [ToL[l3l|. This 
has the advantage that the coupling constant g disappears from the coupled system of equations, as long as gluonic 
vertices are ignored. However, since the diagram containing the three-gluon vertex has a prefactor g 2 but only a 
single ghost propagator in the integral, in Eq. (HO)) there remains an explicit factor g. We will take the coupling at 
the renormalization point /i = [13| g r = 3.5. 

Figure [2J shows the values of the two integrals I\ and Ii separately, for equal ghost and gluon momentum. The 
integral 1\ involving three ghost-gluon vertices approaches the IR value evaluated in Eq. (031 . showing a modest 
dependence on the kinematic configuration. At high momenta it drops off logarithmically, due to the anomalous 
dimension of the ghost propagator. The integral li involving the three-gluon vertex vanishes in the deep IR, due to 
the IR divergence of the gluon energy Q(k) in the denominator, see Eqs. (fTT)) and ([20)1 . and drops off in the UV more 
slowly than the ghost-loop term (l/y/\np instead of 1/lnp), since it contains only one ghost propagator. Figure ^ 
shows the total form factor h(p 2 ; p 2 , xp 2 ) . 

More interesting is the kinematic configuration where the ghost legs have momenta of equal magnitude and the 
gluon momentum is varied, see Fig. [4) One observes that the maximum in the mid-momentum regime gets stronger as 
the gluon momentum approaches zero. If the gluon momentum is very small but not zero (dashed line with x = 10~ 4 
in Fig. 2J) , the form factor is almost indistinguishable from the case of vanishing gluon momentum in the IR and in 
the mid-momentum regime; however, at higher momenta the term involving the three-gluon vertex dominates, and 
we observe again the (lnp)^ 1 / 2 behaviour. 

As explained in the footnote of the introduction, the 3-dimensional Euclidean Yang-Mills theory in Landau gauge is 
equivalent to the Hamiltonian approach in Coulomb gauge in 3 + 1 dimensions using a specific vacuum wave functional, 
which is a decent approximation to the true one in the mid-momentum regime. Therefore we can compare our results 
with lattice data from 3-dimensional Landau gauge [20| , see Fig. [5] The vertex there displays a maximum in the 
mid-momentum regime which is somewhat stronger than in our work; furthermore, the form factor from the lattice 
seems to approach zero or even a small negative value in the IR. Qualitatively similar results have been obtained in 
3- and 4-dimensional Landau gauge both in the continuum [l8j and on the lattice (l9j . 
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FIG. 3. Numerical results for the form factor h(p \p ,xp ) of the ghost-gluon vertex with equal ghost and gluon momentum. 
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FIG. 4. Numerical results for the form factor h(xp 2 ;p 2 ,p 2 ) with equal incoming and outgoing ghost. Explanation see text. 



V. CONCLUSIONS 



We have investigated the ghost-gluon vertex of the Hamiltonian approach to Yang-Mills theory in Coulomb gauge. 
The DSE for the ghost-gluon vertex was solved in one-loop truncation using the non-perturbative ghost and gluon 
propagators obtained previously in the variational approach assuming a bare ghost-gluon vertex. The dressing of the 
ghost-gluon vertex was found to increase the vertex between 15% and 25% in the IR and to vanish asymptotically in 
the UV. Since the gluon propagator obtained in the variational approach with a bare ghost-gluon vertex is qualitatively 
similar to that obtained on the lattice (2(| (which contains the full dressing of vertices) we expect that a fully self- 
consistent solution of the coupled DSEs for the propagators and vertices yield similar results for the ghost-gluon vertex 
as obtained in the present paper. Finally it will be interesting to extend the present studies to finite temperatures (3l| 
to see whether there is a similar change of the dressing of the ghost-gluon vertex at high temperatures as observed in 
the FRG approach in Landau gauge [21 1 . 
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FIG. 5. Comparison of our result for the form factor of the ghost-gluon vertex to lattice data from Ref. [20j at the symmetric 
point, h(p 2 ;p 2 ,p 2 ). 
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